#ifndef __DUNE_COMPOSE__
#define __DUNE_COMPOSE__
// C++ includes
#include<math.h>
// #include <iostream>
// dune-common includes
// #include <dune/common/parallel/mpihelper.hh> // An initializer of MPI
// #include <dune/common/exceptions.hh> // We use exceptions
// #include <dune/common/parametertreeparser.hh>
// dune-grid includes
#include <dune/grid/yaspgrid.hh>
#include <dune/grid/io/file/vtk/vtkwriter.hh>
#include <dune/grid/io/file/vtk/subsamplingvtkwriter.hh>
#include<dune/grid/io/file/vtk/vtksequencewriter.hh>
#include <dune/grid/io/file/gmshreader.hh>
#include <dune/grid/utility/structuredgridfactory.hh>
#if HAVE_UG
#include<dune/grid/uggrid.hh>
#endif
#if HAVE_DUNE_ALUGRID
#include<dune/alugrid/grid.hh>
#include<dune/alugrid/dgf.hh>
#include<dune/grid/io/file/dgfparser/dgfparser.hh>
#endif
// dune-pdelab includes
#include<dune/pdelab/common/function.hh>
#include<dune/pdelab/common/vtkexport.hh>
#include<dune/pdelab/common/quadraturerules.hh>
#include<dune/pdelab/finiteelement/localbasiscache.hh>
#include<dune/pdelab/finiteelementmap/pkfem.hh>
#include<dune/pdelab/finiteelementmap/qkfem.hh>
#include<dune/pdelab/finiteelementmap/qkdg.hh>
#include<dune/pdelab/constraints/common/constraints.hh>
#include<dune/pdelab/constraints/common/constraintsparameters.hh>
#include<dune/pdelab/constraints/conforming.hh>
#include<dune/pdelab/function/callableadapter.hh>
#include<dune/pdelab/gridfunctionspace/gridfunctionspace.hh>
#include<dune/pdelab/gridfunctionspace/gridfunctionspaceutilities.hh>
#include<dune/pdelab/gridfunctionspace/interpolate.hh>
#include<dune/pdelab/gridfunctionspace/vtk.hh>
#include<dune/pdelab/gridoperator/gridoperator.hh>
#include<dune/pdelab/gridoperator/onestep.hh>
#include<dune/pdelab/localoperator/defaultimp.hh>
#include<dune/pdelab/localoperator/pattern.hh>
#include<dune/pdelab/localoperator/flags.hh>
#include<dune/pdelab/localoperator/variablefactories.hh>
#include<dune/pdelab/backend/istl.hh>
#include<dune/pdelab/stationary/linearproblem.hh>
#include<dune/pdelab/instationary/onestep.hh>
#include<dune/pdelab/solver/newton.hh>
#include<dune/pdelab/function/minus.hh>

// dune-[our module] includes
#include <dune/constant-unit/constant.hh>
#include <dune/constant-unit/unit.hh>
#include <dune/equation/equation.hh>
#include <dune/initial/initial.hh>
#include <dune/boundary/boundary.hh>
#include <dune/flux/fluxdeduce.hpp>
#include <dune/compose/adapter/adapter.hh>
#include <dune/compose/wrapper.hpp>
#include <dune/compose/getdt.hpp>
#include <dune/compose/dghyperbolicoperator.hpp>
#include <dune/compose/driver.hpp>
#include <dune/compose/assign-grid.hpp>
#include <dune/compose/assign-equation.hpp>

#endif  //  __DUNE_COMPOSE__
